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1.  Introduction 


The  motivation  for  this  study  is  to  better  understand  the  high-angle-of-attack  flight 
behavior  of  spin- stabilized  projectiles.  Gyroscopically  stabilized  projectiles 
launched  at  steep  angles  for  indirect  fire  experience  large  yaw  of  repose  in  the 
subsonic  regime.  The  uncertainties  associated  with  the  aerodynamics  at  these  high 
angles  of  attack  increase  the  delivery  error  due  to  poor  fire-control  solutions  (i.e., 
understanding  the  relationship  between  the  gun  pointing  angle  and  the  impact 
location).  Guidance  kits  emplaced  on  these  spin- stabilized  vehicles  can  also  induce 
significant  angle  of  attack  to  maneuver.  Accurate  characterization  of  these 
aerodynamics  is  critical  to  removing  miss  distance,  especially  since  control 
authority  is  often  limited  in  these  applications. 

Guns  are  rifled  to  impart  spin  to  projectiles,  which  can  improve  delivery  accuracy 
and  lower  velocity  decay  over  aerodynamically  stabilized  projectiles.  The 
overwhelming  majority  of  spin-stabilized  projectiles  undergo  motion  solely  due  to 
the  external  forces  and  moments  of  gravity  and  rigid  aerodynamic  surfaces  (i.e., 
flight  is  not  controlled  by  actuating  moveable  aerodynamic  surfaces).  The  angle  of 
attack  in  these  applications  is  usually  below  10°.  For  this  reason,  the  literature  on 
spin- stabilized  projectile  aerodynamics  has  focused  on  lower  angles  of  attack.  High 
angle  of  attack  will  be  defined  for  the  purposes  of  this  study  to  be  over  10°.  While 
static  aerodynamics  (e.g.,  drag  force,  lift  force,  and  overturning  moment)  are 
important  to  spin-stabilized  flight,  the  rate-dependent  aerodynamics  (e.g.,  pitch 
damping  moment,  Magnus  moment,  and  roll  damping  moment)  are  also  critical  due 
to  the  role  these  phenomena  play  in  ballistic  flight  stability.1  Wind  tunnel 
techniques  have  been  employed  using  spinning  models2-5  to  directly  measure  forces 
and  moments  but  suffer  from  typical  sting  and  wall  interference  effects  and 
encounter  difficulties  determining  dynamic  aerodynamic  coefficients.  The  spark 
range,6-8  onboard  sensor,9-13  and  high-speed  video14  techniques  have  no  sting  or 
wall  effects  since  the  vehicle  is  in  free  flight.  These  techniques  provide  good 
dynamic  aerodynamics  but  do  not  directly  measure  forces  and  moments  and  feature 
limited  control  over  the  angle  of  attack  range.  Computational  fluid  dynamics 
techniques15  20  have  enjoyed  increasing  success  when  applied  to  spin- stabilized 
projectile  aerodynamics  over  the  past  35  years. 

The  contributions  of  the  present  work  include  defining  the  state  of  the  art  in 
experimental  flight  techniques  and  improving  the  understanding  of  nonlinearities 
in  the  aerodynamics  of  spin- stabilized  projectiles.  To  accomplish  this,  nonlinear 
and  linear  flight  models  were  formulated  with  the  appropriate  aerodynamic  models. 
Experiments  were  designed  and  executed  to  capture  accurate  data  on  models  flying 
at  high  subsonic  to  low  transonic  speeds  with  varying  degrees  of  angle  of  attack 
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from  very  low  (~0°)  to  very  high  (-40°)  using  repeatable  yaw  induction 
techniques.21  Parameter  estimation  algorithms  were  applied  that  resulted  in 
residuals  on  the  order  of  the  measurement  errors  and  subsequently  provided  the 
angle-of-attack  dependency  in  the  aerodynamic  coefficients.  These  experiments 
validated  the  modeling  and  data  analysis  approach  and  ultimately  provide  insight 
into  how  to  improve  delivery  accuracy  of  spin-stabilized  projectiles  flying  at  high 
angle  of  attack. 

2.  Flight  Modeling 


The  flight  modeling  approach  in  this  study  seeks  the  simplest  form  that  represents 
the  relevant  phenomenon.  The  following  equations  govern  the  flight  motion  of 
gyroscopically  stabilized  vehicles  and  can  be  found  in  the  literature.1-22  The  rigid- 
body  kinematics  (Eqs.  1  and  2)  and  dynamics  (Eqs.  3  and  4)  are  obtained 
straightforwardly  from  considering  first  principles.  The  rotational  dynamics 
expressions  used  in  this  study  (Eq.  4)  are  somewhat  complicated  by  the  inclusion 
of  asymmetric  mass  configurations  in  experiments. 
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A  more  active  area  of  research  in  flight  modeling  is  the  aerodynamic  model 
(Eqs.  5-10),  which  is  often  configuration-dependent.  The  aerodynamic  model  used 
in  this  work  includes  nonlinearities  in  the  static  forces  and  moments,  damping 
moments,  and  nonlinear  Magnus  (side)  moments.  Aerodynamic  coefficients 
depend  on  Mach  number  and  total  angle  of  attack.  Variation  with  aerodynamic  roll 
angle  is  neglected  (e.g.,  assume  Cn  =  Cma)  due  to  the  configurational  symmetry 

in  the  spin- stabilized  vehicle  shape.  While  polynomials  are  used  for  most 
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aerodynamic  phenomena  a  table  is  implemented  to  capture  the  large  nonlinearity  in 


the  Magnus  moment  with  angle  of  attack. 

X  =  -Qs(cXo  +  Cx_a2S\n2a).  (5) 

Y  =  -QS  (cNas\n  p  +  CNa3s\r\3  p  -  CYpa^s\r\a).  (6) 

Z  =  -QS  (cWaSina  +  CNa3  sin3  a  +  CYpJ^s\r\  p}.  (7) 

L  =  QSDClp (8) 
M  =  QSD  (cma  sin  a  +  Cma,  sin3  a  +  Cmq  +  CnJ/?)  f^sin  /?).  (9) 

N  =  QSD  (-Cma  sin /?  -  Cma 3  sin3/?  +  Cm<j g  +  CnJa)  f^sin a).  (10) 


Multiple  spark  range  flights  in  different  angle  of  attack  ranges  yields  the  best 
information  concerning  the  nonlinear  aerodynamics.  Various  masses  offset  from 
the  body  axis  of  symmetry  were  placed  within  the  projectiles  to  repeatably  induce 
several  magnitudes  of  yaw.  Linearization  of  the  flight  motion1,22  provides  analytical 
relationships  for  the  vehicle  angular  motion.  This  angular  motion  has  components 
due  to  nutation,  precession,  and  trims.  Trims  produce  angular  oscillations  at  the 
spin  rate.  Spin  rate  is  given  by  gun  rifling,  launch  velocity,  axial  moment  of  inertia, 
and  roll  damping  moment  for  this  application.  The  frequencies  of  nutation  and 
precession  can  be  obtained  from  the  aerodynamics,  mass  properties,  flight  states, 
and  atmospheric  conditions,  as  given  in  Eq.  11. 
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Equations  relating  the  tilt  angle  (angle  between  the  principal  inertial  axis  and  the 
body  axis  of  symmetry)  and  the  maximum  total  angle  of  attack23  can  be  exploited 
for  experiment  design.  Equation  12  was  used  to  size  different  mass  inserts  for  the 
projectiles  during  firings.  The  magnitude  of  the  maximum  total  angle  of  attack 
increases  significantly  for  high  spin  rates,  which  are  usually  encountered  for 
gyroscopically  stabilized  applications. 

■w  -  <12> 

L  COft-COp  J 


3.  Spark  Range  Firings 

Fifteen  vehicles  were  flown  through  the  spark  range  with  7  different  mass  insert 
configurations.  The  mass  properties  were  obtained  using  a  special  procedure 
crafted  to  determine  the  full  moment  of  inertia  tensor  on  precision  measurement 
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devices.  Mass  properties  of  the  vehicles  are  provided  in  Table  1.  The  tilt  axis  is 
related  to  the  moments  of  inertia  ( IXY  =  Iyy  ^'xx  tan  2y).  The  external  shape  of  the 
projectile  was  characteristic  of  an  artillery  application. 

Table  1  Mass  properties  of  flight  vehicles 


No.  of 
Flights 

Mass 

(kg) 

Longitudinal 
center  of 
gravity3 
(mm) 

Lateral 
center  of 
gravity3 
(mm) 

Axial 

inertia 

(kg-m2) 

Trans. 

inertia 

(kg-m2) 

Cross¬ 

axis 

inertia 

(kg-m2) 

2 

25.3921 

0.559499 

0 

0.114429 

1.166983 

0 

1 

26.8250 

0.571576 

0.001016 

0.116729 

1.311924 

0.001835 

4 

27.3153 

0.575840 

0.000787 

0.117523 

1.311924 

0.007454 

1 

27.9322 

0.578321 

0.001886 

0.119091 

1.311924 

0.008985 

1 

28.2588 

0.580697 

0.002273 

0.119250 

1.311924 

0.011786 

4 

28.8394 

0.582118 

0.003000 

0.120505 

1.318907 

0.012992 

2 

29.2658 

0.587317 

0.004114 

0.120707 

1.370948 

0.019349 

a  From  nose 


A  155-mm-diameter  rifled  gun  launched  the  vehicles.  A  half  muzzle  break  was  not 
used  to  induce  yaw  in  an  attempt  to  isolate  angular  motion  due  to  asymmetric  mass. 
The  gun  was  mounted  in  a  tracked  vehicle  and  positioned  with  the  muzzle  about 
2  m  from  the  spark  range  entrance.  The  propelling  charge  used  1.18  kg  of  Ml 
propellant  with  about  30  gr  of  black  powder  for  the  primer.  These  conditions 
resulted  in  launch  speeds  in  the  low  transonic  to  high  subsonic  regime. 

A  series  of  preliminary  firings  were  conducted  with  high-speed  photography  (fixed 
and  flight  follower  using  a  rotating  mirror),  radar,  yaw  cards,  and  gun  pressure 
gauges  to  establish  the  charge,  verify  structural  integrity,  assess  flight  stability,  and 
validate  the  flight  path  prior  to  spark  range  firings.  These  firings  were  especially 
critical  since  the  vehicles  must  fly  through  an  approximately  2-  x  2-m  window  at 
each  station  in  the  spark  range.  The  gravity  drop  is  significant  at  this  Mach  number, 
and  the  swerving  motion  is  large  for  high-angle-of- attack  experiments.  A 
representative  flight  follower  dataset  on  a  configuration  that  reached  about  25° 
angle  of  attack  is  provided  in  Fig.  1. 
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Fig.  1  Flight  follower  of  spin-stabilized  vehicle  flying  to  approximately  25°  angle  of  attack 

The  spark  range  technique  has  been  applied  with  success  over  many  years  to  study 
the  flight  behavior  of  different  vehicle  types  such  as  aerodynamically  stabilized 
projectiles  (e.g.,  tank  ammunition24  and  mortars),  aircraft,6  re-entry  vehicles,25 
guided  munitions,8'26  as  well  as  spin-stabilized  projectiles.7  The  spark  range  at  the 
US  Army  Research  Laboratory’s  Transonic  Experimental  Facility27  features  25 
measurement  stations  in  an  enclosed  building  more  than  250  m  long.  Each  station 
is  equipped  with  a  light  screen,  delay  unit,  spark  source,  surveyed  fiducials,  and  2 
orthogonal  cameras  for  capturing  shadowgraphs.  As  the  vehicle  flies  past  each 
station,  the  light  screen  triggers,  which  sets  off  the  high-intensity,  short-duration 
spark  source  after  the  appropriate  delay  to  obtain  the  orthogonal  shadowgraphs. 
Postprocessing  of  the  shadowgraphs  and  trigger  times  yields  the  center-of-gravity 
location  and  Euler  angles,  which  along  with  meteorological  data  are  used  to  obtain 
aerodynamic  coefficients  in  a  parameter  estimation  algorithm. 

4.  Parameter  Estimation  Algorithm 

The  maximum  likelihood  method2528’29  was  implemented  to  determine 
aerodynamics  from  these  spark  range  firings.  This  is  a  nonlinear  least-squares 
technique  that  seeks  to  find  parameters  that  minimize  a  logarithmic  function 
containing  the  residual  between  the  experimental  measurements  and  the  model 
calculations.  The  parameter  vector  for  this  problem  included  the  aerodynamic 
coefficients,  initial  conditions,  and  cross-axis  moment  of  inertia.  Magnus  force  was 
not  included  in  the  parameter  estimation.  The  Magnus  moment  table  was  adjusted 
using  linear  and  cubic  terms  (Eq.  13). 
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Initial  guesses  for  these  parameters  were  used  to  integrate  both  the  nonlinear  flight 

d  x 

model  (Eqs.  1-10)  and  differential  equations  for  a  Jacobian  (— )  to  the  end  of  the 

time  series.  A  Newton-Raphson  routine  was  applied  to  optimize  the  likelihood 
function  at  times  containing  experimental  data.  Equation  14  was  used  in  this 
method  to  update  the  values  in  the  parameter  vector. 


A  Cl  —  ( YnM  dxc,iT  ira-1  d*c,i  ,  rj\  dxc/ 

-w  +  V)  U= iii- 


(14) 


This  procedure  was  iterated  with  the  adjusted  parameters  until  convergence  was 
achieved. 


5.  Results 


Representative  shadowgraphs  are  provided  in  Fig.  2.  These  data  were  collected  at 
the  same  instant  in  time  in  the  horizontal  (Fig.  2,  bottom)  and  vertical  (Fig.  2,  top) 
planes.  In  each  image  the  blurred  shape  is  the  actual  vehicle,  and  the  sharp  shape  is 
the  shadow  cast  onto  the  panel  by  the  spark  source.  At  this  station  the  vehicle  is 
pitched  up  about  40°  with  a  yaw  angle  near  0°.  Note  the  pin  inserted  into  the  base 
to  measure  roll  angle  along  with  some  wake  flow  structures  trailing  the  body. 
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Fig.  2  Orthogonal  shadowgraphs 
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Figures  3-6  show  the  measured  and  reconstructed  motion  for  one  of  the  flights.  In 
these  plots,  the  dots  represent  the  experimental  measurements  and  the  curves 
represent  the  model  calculations  yielded  by  the  parameter  estimation  algorithm. 
The  horizontal  and  vertical  center-of-gravity  locations  are  illustrated  in  Fig.  3.  The 
action  of  gravity  is  apparent  in  the  data  in  the  vertical  plane.  The  horizontal  data 
increase  significantly  with  downrange  location  due  to  jump.  This  high-amplitude 
swerving  motion  is  important  because  the  accuracy  in  the  estimated  aerodynamic 
normal  force  increases  with  magnitude  of  the  swerve.  A  swerving  period 
(approximately  2  cycles  in  200  m  of  downrange  travel)  is  also  evident  in  the 
horizontal  data. 


Fig.  3  Center-of-gravity  motion 


The  rolling  motion  is  captured  in  Fig.  4.  These  data  are  strongly  linear  with  distance 
travelled  because  the  launch  spin  rate  is  around  90  Hz  and  roll  decelerates  slightly 
due  to  aerodynamic  damping.  This  vehicle  rolls  about  70  times  over  200  m  of 
downrange  travel.  For  some  of  the  high-angle-of-attack  flights  a  significant  amount 
of  roll  data  were  missing  due  to  occlusions  (i.e.,  the  base  prevented  the  pin  from 
being  in  view),  which  limited  the  roll  analysis.  This  was  not  problematic  for  the 
overall  flight  analysis  because  the  roll  motion  is  predictable  for  spin-stabilized 
projectiles  flying  these  distances. 
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The  pitch  and  yaw  motion  for  this  flight  is  provided  in  Figs.  5  and  6.  Figure  5 
illustrates  the  tri-cyclic  motion  undergone  by  this  vehicle.  The  nutation  mode  acts 
at  the  lowest  frequency  (~2  Hz),  followed  by  the  precession  mode  at  a  slightly 
higher  frequency  (~6  Hz),  and  finally  the  trim  mode,  which  is  at  the  spin  rate 
(~90  Hz).  Each  of  these  modes  is  clearly  visible  by  separate  frequencies  in  the  data. 
The  data  in  Fig.  5  is  plotted  differently  in  Fig.  6  to  show  the  complexity  of  this 
motion  when  viewed  from  the  fixed-plane  body  reference  and  because  this  vehicle 
flies  from  about  0°  to  15°  total  angle  of  attack.  Collecting  data  over  angle-of-attack 
ranges  such  as  these  is  critical  to  investigating  the  nonlinear  aerodynamics. 

20r 
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Fig.  6  Fixed-plane  angular  motion 

Four  firings  of  identical  configurations  with  different  initial  roll  orientations  were 
conducted  and  data  analysis  performed  to  identify  the  relationship  between  the 
internal  mass  orientation  and  the  swerve  motion  direction.  These  results  are 
presented  in  Fig.  7.  Each  flight  corresponds  with  a  particular  swerve  direction  (up, 
right,  down,  and  left)  as  shown  in  the  legend.  The  horizontal  data  given  on  the  left 
of  the  figure  is  viewed  from  above  looking  toward  the  ground.  The  left  flight 
horizontal  position  increases  with  downrange  distance,  and  the  right  flight 
horizontal  location  decreases  with  downrange  distance.  The  up  and  down  flights 
remain  approximately  on  the  line  of  fire.  The  total  span  of  swerve  motion  is 
approximately  1.5  m  depending  on  the  initial  roll  orientation.  The  oscillation  in  the 
swerve  due  to  nutation  and  precession  angular  motion  is  apparent  in  the  horizontal 
data.  Inspection  of  the  vertical  data  at  the  right  of  the  figure  shows  a  higher 
trajectory  for  the  up  flight  and  a  lower  trajectory  for  the  down  flight.  The  left  flight 
is  slightly  higher  than  the  right  flight  due  to  a  higher  launch  elevation  angle.  The 
angle  of  attack  histories  for  these  flights  with  identical  configurations  was  within 
about  20%,  which  indicates  that  this  technique  of  using  internal  mass  asymmetries 
for  yaw  induction  is  fairly  repeatable. 
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Fig.  7  Relationship  between  internal  mass  orientation  and  swerve  direction 


To  obtain  these  swerve  directions  for  spin- stabilized  projectiles,  the  additional  mass 
near  the  base  must  be  90°  counter-clockwise  (when  viewed  from  the  base)  at  launch 
from  the  desired  swerve  direction.  As  an  example,  if  the  desired  swerve  direction 
is  up  (i.e.,  12:00)  then  at  launch  the  mass  insert  near  the  base  must  be  oriented  to 
the  left  (i.e.,  9:00)  when  viewed  from  behind  the  gun.  Gun  tube  length,  rifling,  and 
load  roll  orientation  must  be  known  to  prescribe  the  swerve  direction  in  this 
manner. 

The  summary  flight  motion  from  all  firings  is  provided  in  Fig.  8.  The  data  were 
translated  to  better  visualize  the  horizontal  and  vertical  center  of  gravity.  The  spread 
in  the  horizontal  and  vertical  center-of-gravity  motion  is  about  2  m  at  200  m 
downrange,  which  is  about  the  size  of  the  spark  range  instrumentation  windows. 
This  suggests  that  experiments  were  designed  to  maximize  the  observation  of  flight 
behaviors  subject  to  experimental  limitations.  The  Euler  angles  qualitatively  show 
the  relationship  between  angular  motion  amplitude  and  the  tilt  axis. 
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Fig.  8  Summary  flight  motion 

Overall,  the  model  matches  the  measurements  to  within  the  experimental  error, 
which  validates  the  aerodynamic  model  (Eqs.  5-10),  kinematic  and  dynamic 
equations  (Eqs.  1-4),  and  parameter  estimation  algorithm.  There  is  larger 
uncertainty  in  the  results  at  higher  angle  of  attack.  Measurement  errors  may  be 
larger  at  higher  angle  of  attack  due  to  perspective  changes  of  the  camera/spark 
source  and  vehicle/shadowgraph.  In  addition,  modeling  errors  may  be  larger  due  to 
neglecting  nonlinear  variation  with  angle  of  attack  in  the  roll  damping  moment  and 
pitch  damping  moment. 

The  total  angle  of  attack  histories  for  all  flights  are  shown  in  Fig.  9.  Experiments 
successfully  captured  a  wide  spectrum  (0°-40°)  of  total  angle  of  attack  that  is 
critical  to  examining  nonlinearities  in  aerodynamics  with  angle  of  attack.  Most 
flights  appears  stable;  however,  2  of  the  flights  at  the  largest  angle  of  attack  indicate 
growth  in  the  angular  motion  amplitude  with  downrange  distance.  The  oscillation 
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frequency  in  the  angular  motion  is  dictated  primarily  by  the  static  pitching  moment. 
The  slight  shifting  in  the  frequencies  of  the  angular  motion  illustrated  in  Fig.  9  is 
likely  attributable  to  spread  in  the  Mach  histories  of  the  flights  in  addition  to 
potential  variation  in  the  static  pitching  moment  with  angle  of  attack  and  Mach 
number.  The  highest  frequency  in  this  data  is  at  the  spin  rate  due  to  the  offset  mass 
trim.  The  2  lowest  angle  of  attack  flights  do  not  exhibit  this  behavior  due  to  the 
absence  of  an  asymmetric  mass  for  that  configuration  (see  Table  1). 


Fig.  9  Summary  angle  of  attack  histories 


The  relationship  between  the  maximum  total  angle  of  attack  and  the  tilt  axis 
obtained  from  this  flight  analysis  and  theory23  (Eq.  12)  is  given  in  Fig.  10.  The 
linear  term  in  a  least-squares  fit  to  the  experimental  data  is  also  shown.  These 
results  extend  those  shown  by  Rollstin  and  Hodapp12  and  imply  that  the  analytical 
equation23  is  valid  for  designing  experiments  such  as  those  conducted  in  this  study. 
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Fig.  10  Relationship  between  maximum  total  angle  of  attack  and  tilt  axis 

The  aerodynamics  that  result  from  the  flight  analysis  are  presented  next.  The 
average  Mach  number  for  all  flights  was  0.77  with  a  standard  deviation  of  0.03. 
The  moment  reference  location  is  3.755  cal  from  the  nose. 

The  axial  force  coefficient  variation  with  angle  of  attack  for  individual  flights  and 
a  least- squares  fit  to  all  the  data  is  provided  in  Fig.  11.  The  zero- yaw  axial  force 
obtained  in  this  study  compares  favorably  with  free-flight  and  wind  tunnel  results 
on  similar  spin- stabilized  shapes  and  Mach  numbers  in  Whyte7  and  Kahn  et  al.10 
Axial  force  increases  approximately  70%  from  0°  to  25°  total  angle  of  attack.  There 
is  spread  in  the  results  at  a  particular  aerodynamic  angle.  At  small  aerodynamic 
angles  this  scatter  can  be  attributed  to  typical  physical  variations  associated  with 
manufacturing  tolerances  and  the  launch  event. 
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Fig.  11  Axial  force  coefficient 

Spread  in  the  data  increases  with  aerodynamic  angle,  which  indicates  that  a 
continuous  polynomial  does  not  represent  the  underlying  flow  phenomena  at  high 
angle  of  attack.  There  is  little  data  in  the  literature  on  the  axial  force  for  spin- 
stabilized  projectiles  at  high  angle  of  attack  on  which  to  base  a  modeling  approach. 
Flight  analysis  was  conducted  using  a  table  (as  done  for  the  Magnus  moment  in  this 
study)  rather  than  a  polynomial  for  the  axial  force  aerodynamic  model,  but  the 
results  did  not  improve. 

Obtaining  accurate  axial  force  at  high  angle  of  attack  is  challenging  in  these 
experiments  for  2  primary  reasons.  First,  aerodynamic  forces  and  moments  affect 
the  accelerations  of  the  vehicle,  but  the  spark  range  measurements  used  for 
parameter  estimation  are  positions  (e.g.,  center-of-gravity  location  and  Euler 
angles).  Nonlinearities  in  the  aerodynamics  with  angle  of  attack  may  be  low-pass 
filtered  as  a  result  of  integration  from  acceleration  (i.e.,  aerodynamic  force  and 
moment)  to  position. 

The  second  difficulty  is  related  to  the  measurements  for  vehicles  flying  at  high 
angle  of  attack.  For  typical  spark  range  flights  at  low  angle  of  attack,  matching  the 
time-stamped  downrange  center-of-gravity  location  provides  the  zero-yaw  axial 
force.  Matching  the  pitch  and  yaw  angles  permits  assessment  of  the  yaw  component 
of  axial  force,  which  at  low  angle  of  attack  is  often  small  compared  with  the  zero- 
yaw  axial  force.  Matching  the  small  swerve  motion  (primarily  the  horizontal  and 
vertical  center-of-gravity  location)  enables  some  evaluation  of  the  normal  force.  In 
this  situation  the  flight  resembles  a  point-mass  model.  At  high  angle  of  attack  the 
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6-degree-of-freedom  model  is  much  more  in  play.  Components  of  axial  force  as 
well  as  normal  force  contribute  to  the  time- stamped  downrange  center-of-gravity 
location.  Nonlinearities  with  angle  of  attack  in  the  axial  force  and  normal  force  are 
more  important  to  matching  the  center-of-gravity  measurements,  so  properly 
estimating  the  aerodynamic  moments  to  match  angular  motion  is  extra  critical.  This 
situation  increases  the  coupling  between  translational  and  rotational  dynamics  of 
the  vehicle.  In  practice,  for  some  of  the  highest  angle-of-attack  flights  when  the 
angular  motion  was  well-fit  by  a  set  of  aerodynamic  moments,  the  axial  force  and 
normal  force  could  be  traded  for  an  equivalent  fit  to  the  center-of-gravity  motion. 
Moreover,  the  highest  angle-of-attack  flights  were  unstable,  which  was  reflected  in 
the  numerical  stability  of  the  flight  model  and  complicated  the  analysis. 

Flying  multiple  identically  shaped  vehicles  at  a  given  Mach  number  over  numerous 
angle-of-attack  (including  high)  ranges  arguably  provides  a  “best”  set  of  spark 
range  data  for  obtaining  aerodynamic  normal  force  (i.e.,  high  signal-to-noise  ratio 
in  lateral  center-of-gravity  motion  data).  These  data  are  plotted  in  Fig.  12.  The 
results  suggest  a  mild  (positive  cubic  term)  nonlinearity  with  angle  of  attack  up  to 
40°.  Wind  tunnel  experiments  with  spinning  models  of  similar  geometry  around 
Mach  0.3  also  showed  primarily  linear  behavior  (small  positive  cubic  term)  in  the 
normal  force  up  to  70°,  which  were  also  largely  independent  of  nondimensional 
spin  rate. 1  The  data  in  Fig.  12  also  agree  with  relevant  free-flight  results  by  Whyte.7 


Fig.  12  Static  normal  force  coefficient 
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Figure  13  gives  the  static  pitching  moment  coefficients.  Similar  to  the  normal  force 
results,  the  static  pitching  moment  is  relatively  linear  (small  negative  cubic  term) 
to  40°.  These  data  agree  with  related  experimental  results  presented  by  McCoy,3 
Whyte,7  and  Kahn  et  al.10  McCoy3  indicates  an  increased  nonlinearity  with  angle 
of  attack  around  30°. 


Figure  14  provides  the  pitch  damping  moment  (left)  and  roll  damping  moment 
(right)  results  as  a  function  of  the  average  Mach  number  of  each  flight  since 
variation  with  angle  of  attack  was  neglected.  The  pitch  damping  moment  features 
some  scatter  typical  of  spin-stabilized  projectiles  due  to  the  ambiguity  between 
pitch  damping  moment  and  Magnus  moment  from  spark  range  observations 
explained  elsewhere.1'30  Pitch  damping  moments  compare  favorably  with  Whyte7 
and  Kahn  et  al. 10 

Roll  damping  moment  was  obtained  with  some  ease  because  the  roll  motion  is 
relatively  decoupled.  There  was  some  uncertainty,  as  described,  due  to  occluded 
roll  measurements,  especially  at  high  angle  of  attack.  The  roll  damping  moment 
results  are  repeatable  and  support  those  found  by  Whyte.7 
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Fig.  14  Pitch  and  roll  damping  moment  coefficients 


Magnus  moment  coefficients  at  an  average  Mach  number  of  0.77  are  provided  in 
Fig.  15.  The  nutation  (negative  value)  and  precession  (positive  value)  stability 
bounds  determined  from  analytical  expressions  based  on  linearizing  the  equations 
of  motion1,30  were  added  to  the  plot  as  the  dashed  lines.  The  Magnus  moment  is 
highly  nonlinear  with  angle  of  attack.  At  low  angles  a  nutation  instability  is 
predicted.  As  the  angle  of  attack  grows  due  to  this  instability,  the  Magnus  moment 
remains  in  a  stable  region  between  about  3°  and  6°.  Above  6°  the  Magnus  moment 
hugs  the  precession  stability  bound.  These  results  suggest  that  this  vehicle  has  a 
weak  precession  instability  at  moderate  and  high  angle  of  attack  so  that  flight 
perturbations  that  produce  total  angle  of  attack  above  6°  could  result  in  constant  or 
increased  total  angle  of  attack.  The  notion  of  linear  stability  is  challenged  at  these 
high  angles  of  attack  because  the  solution  assumes  low  angle  of  attack.  Better 
representations  of  stability  may  provide  some  insight  into  how  stability  bounds 
change  with  angle  of  attack  or  how  quickly  angular  motion  changes  depending  on 
distance  from  the  stability  bound. 

As  discussed  earlier,  there  are  challenges  in  characterizing  angle-of-attack 
nonlinearities  in  aerodynamics  from  free-flight  experiments.  There  is  uncertainty 
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in  the  Magnus  moment  at  low  angles  because  it  is  difficult  to  fly  vehicles  at  very 
low  angle  of  attack  due  to  launch  disturbances,  internal  and  external  vehicle 
asymmetries,  and  atmospheric  effects  such  as  wind.  Combining  this  limited  low 
angle  of  attack  free-flight  dataset  and  the  current  modeling  approach  yields  a 
Magnus-induced  nutation  instability  at  small  angles. 
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Fig.  15  Magnus  moment  coefficient 

The  increase  in  Magnus  moment  with  angle  of  attack  is  supported  in  the  literature 
on  spinning  bodies  of  revolution  at  subsonic  Mach  number.3,5’7,10  Whyte5  indicated 
a  negative  Magnus  moment  at  low  angles.  McCoy3  and  Kahn  et  al. 10  were  the  only 
studies  above  8°  angle  of  attack.  The  nonlinearities  in  Magnus  moment  with  angle 
of  attack  do  not  agree  between  these  studies  but  that  could  be  due  to  differences  in 
Mach  number  or  external  shape.  McCoy3  observed  large  out-of-plane  oscillations 
in  spinning  wind  tunnel  models  for  angle  of  attack  between  30°  and  150°.  This 
phenomenon  required  averaging  of  the  noisy  gauge  balance  data  and  increased  the 
uncertainty  in  the  subsequent  aerodynamic  coefficients  such  as  Magnus  moments 
in  this  angle  of  attack  range.  Unsteady  vortex  shedding,  similar  to  that  found  on 
slender  bodies  at  high  angle  of  attack31  or  long  circular  cylinders  in  a  cross-flow 
(i.e.,  Karman  vortex  street),  was  offered  to  explain  this  behavior.  Unsteady  effects 
such  as  these  were  investigated  by  DeSpirito32  and  need  further  study  to  close  gaps 
in  understanding  the  interaction  between  the  flow  and  the  vehicle  motion.  Time- 
dependencies  can  be  incorporated  into  the  aerodynamic  model  (Eqs.  5-10)33  at  the 
appropriate  aerodynamic  angles  to  represent  these  phenomena. 
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6.  Conclusions 


This  study  was  the  first  to  characterize  the  aerodynamics  of  spin-stabilized  vehicles 
in  this  Mach  regime  at  high  angle  of  attack  from  free-flight  spark  range  data.  The 
key  findings  and  contributions  to  the  state  of  the  art  are  separated  into  3  areas: 
experimental  techniques,  theory  validation,  and  understanding  of  spin-stabilized 
vehicle  flight  behavior  at  high  angle  of  attack. 

Experiments  were  successful  in  demonstrating  an  approach  for  obtaining  the 
desired  high  angle-of-attack  flight  behaviors  while  staying  within  spark  range 
instrumentation  windows.  Vehicle  physical  characteristic  modeling  and 
measurements,  preflight  simulations,  and  preliminary  firings  with  high-speed 
photography,  radar,  yaw  cards,  and  gun  pressure  gauges  were  essential  to  designing 
the  experiments  and  verifying  the  launch  and  flight  system.  Flights  were  repeatable, 
which  enabled  experiments  to  be  performed  up  to  40°  angle  of  attack.  The 
consistency  of  results  allowed  the  relationship  between  launch  configuration  roll 
orientation  and  swerve  direction  to  be  defined.  Accurate  raw  data  (e.g.,  asymmetric 
mass  properties  and  shadowgraphs)  were  collected.  Measurement  errors  due  to  the 
effect  of  perspective  errors  in  the  shadowgraph  photography  system  at  high  angle 
of  attack  should  be  studied. 

Theoretical  expressions  relating  the  vehicle  tilt  axis  to  the  maximum  total  angle  of 
attack  were  validated  up  to  40°.  This  lends  credibility  to  applying  these  expressions 
for  experimental  design  or  performance  analysis.  Nonlinear  stability  theory  should 
be  explored  further  for  spin- stabilized  vehicles  because  linear  stability  bounds  are 
inadequate  at  high  angle  of  attack  and  do  not  provide  suitable  metrics.  Flight 
analysis  mainly  achieved  the  expected  measurement  errors,  which  validates  the 
kinematics  and  dynamics  with  internal  asymmetries  (Eqs.  1-4),  aerodynamic 
modeling  (Eqs.  5-10),  and  parameter  estimation  algorithms.  This  study  revealed 
some  limitations  associated  with  obtaining  aerodynamics  from  spark  range  flights. 
Nonlinearities  in  the  aerodynamics  with  angle  of  attack  are  difficult  to  accurately 
assess  because  the  aerodynamic  forces  and  moments  that  appear  in  the  acceleration 
equations  can  be  low-pass  filtered  during  analysis  to  the  measurements  of  position. 
While  some  of  these  problems  may  be  remedied  with  advanced  parameter 
estimation  algorithms  and  aerodynamic  modeling,  there  are  fundamental 
limitations  associated  with  the  measurement  basis. 

This  study  has  extended  the  understanding  of  spin-stabilized  vehicle  flight  in  the 
high  subsonic  to  low  transonic  Mach  regime  for  high  angle  of  attack.  This 
knowledge  can  be  applied  to  improve  the  accuracy  of  indirect  fires  for  unguided 
(e.g.,  fire  control)  and  guided  (e.g.,  fire  control,  guidance  and  flight  control 
algorithms)  munitions  launched  at  steep  elevation  angles  that  encounter  high  yaw 
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of  repose.  To  improve  prediction  of  the  impact  point  under  these  conditions  it  is 
critical  to  accurately  assess  the  axial  force  and  normal  force.  The  axial  force  dictates 
the  decay  of  velocity  and  subsequent  position  history  along  the  flight  path.  The 
normal  force  drives  the  drift  of  spin-stabilized  projectiles  (to  the  right  from  behind 
the  gun)  off  the  line  of  fire.  This  study  indicates  significant  nonlinearity  in  axial 
force  with  angle  of  attack  that  may  not  be  represented  with  a  polynomial. 
Quantifying  this  nonlinearity  and  the  underlying  aerodynamic  model  structure  are 
critical  to  impact  point  prediction.  The  results  of  this  study  suggest  that  normal 
force  is  relatively  linear  in  angle  of  attack  up  to  40°,  so  devoting  resources  to 
quantify  this  aerodynamic  coefficient  at  high  angle  of  attack  may  not  be  as 
important  as  axial  force. 

The  discussion  of  the  relationship  between  the  Magnus  moment  and  stability 
underscored  how  spin-stabilized  vehicles  fly  within  a  range  of  angle  of  attack.  Extra 
effort  should  be  taken  in  characterizing  the  aerodynamics  across  this  angle-of- 
attack  range.  The  Magnus  moment,  pitch  damping  moment,  and  static  pitching 
moment  primarily  contribute  to  defining  this  angle-of-attack  range.  This  study 
showed  that  for  flat-fire  trajectories  (i.e.,  negligible  yaw  of  repose)  under  nominal 
conditions,  the  angle-of-attack  range  is  between  about  3°  and  6°.  For  steep  elevation 
launch  (-70°)  at  high  charge,  the  upper  end  of  the  angle-of-attack  range  can  be  as 
much  as  35°  due  to  yaw  of  repose  (based  on  low  angle-of-attack  aerodynamics). 
The  static  pitching  moment  was  mainly  linear  and  the  Magnus  moment  was  highly 
nonlinear  with  angle  of  attack  in  this  study.  The  pitch  damping  moment  and  roll 
damping  moment  were  assumed  independent  of  angle  of  attack  in  this  work. 

Overall,  these  results  suggest  that  to  improve  delivery  accuracy  for  spin-stabilized 
projectiles  flying  at  high  angle  of  attack  due  to  steep  elevation  launch,  it  is  of 
primary  importance  to  quantify  a  set  of  highly  nonlinear  aerodynamics  (axial  force 
and  Magnus  moment)  and  of  secondary  importance  to  quantify  a  set  of  largely 
linear  aerodynamics  (normal  force  and  static  pitching  moment)  from  about  3°  to 
40°.  Recommendations  for  the  pitch  damping  and  roll  damping  moments  are 
outside  the  scope  of  the  present  work  but  have  been  investigated  by  DeSpirito.32 
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